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Abstract 

We study dynamics of a coupled system consisting of the 3D Na- 
vier-Stokes equations which is hnearized near a certain Poiseuihe type 
flow in an (unbounded) domain and a classical (possibly nonlinear) 
elastic plate equation for transversal displacement on a flexible flat 
part of the boundary. We first show that this problem generates an 
evolution semigroup St on an appropriate phase space. Then under 
some conditions concerning the underlying (Poiseuille type) flow we 
prove the existence of a compact finite-dimensional global attractor 
for this semigroup and also show that St is an exponentially stable Co- 
semigroup of linear operators in the fully linear case. Since we do not 
assume any kind of mechanical damping in the plate component, this 
means that dissipation of the energy in the fluid flow due to viscosity 
is sufficient to stabilize the system. 

Keywords: Fluid-structure interaction, linearized 3D Navier-Stokes 
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1 Introduction 

Let O C M'^ be a (possibly unbounded) domain with a sufficiently smooth 
boundary dO. We assume that dO = U 5, where n S = 0, 

{x= (xi; X2; 0) : x' = (xi; X2) G M^} 

is bounded in M? and has the smooth contour F = dCl. We refer to Assump- 
tion [2?T] below for further hypotheses concerning the domain O. 
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Let ao(x) = (aKx); aQ{x); Uq^x)) be a smooth bounded field defined on 
O such that divao = 0, (n, ap) = on dO (n is the exterior normal to dO, 
n = (0; 0; 1) on Q) and A = A(x) be a bounded measurable 3x3 matrix, 
X € O. We introduce a linear first order operator Lq of the form 

Lov = {ao,V)v + Av (1) 

and consider the following linear Navier-Stokes equations in O for the fluid 
velocity field v = v{x,t) = {v^{x,t);v'^{x,t);v^{x,t)) and for the pressure 
p{x,t): 

vt - uAv + Lqv + Vp = Gf{t) in C'x(0,+oo), (2) 

divt; = in Ox(0,+oo), (3) 

where > is the dynamical viscosity, Gf{t) is a volume force (which may 
depend on t). We supplement ([2]) and ([3]) with the (non-slip) boundary 
conditions imposed on the velocity field v = v{x,t): 

v = on S; v = {v^;v'^;v^) = {0;0;ut) on f]. (4) 

Here, as in [13], u = u{x,t) is the transversal displacement of the plate 
occupying il. and satisfying the following equation 

uu + A^u + T{u) = Gpi{t) +p\n in nx{0,oo), (5) 

where J-{u) is a nonlinear (feedback) force (see Assumption 14. 1 1 below) . p is 
the pressure from ([2]), Gpi{t) is a given external (non-autonomous) load. We 
refer to [14] and to the references therein for some discussion of this plate 
model and for an explanation of the structure of the force exerted by the 
fluid on the plate. 

We impose clamped boundary conditions on the plate deflection 



du 
on 



= (6) 

an 



and supply ©-([B]) with initial data of the form 

v{0) = vq, u{0) = uq, ut{0) = ui. (7) 

If we assume that the velocity field v decays sufficiently fast as — >• +oo 
and X G O, then ([3|) and (jH) imply the following compatibility condition 

ut{x',t)dx' = forallt>0, (8) 

a 

which can be interpreted as preservation of the volume of the fluid. 

Below (see Definitions 13.11 and 14. 2p we define a solution to ©-([SI) as 
a pair (f ; u) satisfying some variational type relation. If the pair {v; u) is 
already determined, then (at least formally) we can find Vp in O and the 
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trace of p on from and . Thus the pressure p is uniquely defined by 
{v;u). 

The main example which we have in mind is the Poiseuille flow (see, 
e.g., [6] for some details). In this case we deal with the domain 

O = {{xi-X2; X3) : (x2; X3) G ^ c xi £ M}, (9) 

where =^ is a domain in M^, and the Poiseuille velocity field has the form 
ao = (a(x2; X3); 0; 0), where a(x2;x3) solves the elliptic problem 

i/Aa = -A; in ^, a = on d^g, (10) 

where A; is a positive parameter. Linearization of the nonlinear Navier-Stokes 
equations around the flow oq gives us the model with 

Lqv = {aQ,V)v + {v,V)ao. (11) 

There are two important special cases of the choice of ^ in ([9]): (i) ^ is a 
bounded domain in (the Poiseuille flow in a cylindrical tube) and (ii) a 
flow between two parallel planes. In the latter case 

^ = {{x2;x3):x2eR,X3£{-h,0)}, 0(^2; X3) = -^(/i + 2:3). (12) 

Another possibility included in the framework above is the Oseen mod- 
ification of dl]) (see, e.g., [6]). In this case Lqv = Udx^v, where U is the 
parameter which has the sense of the speed of the unperturbed flow mov- 
ing along the xi-axis. This corresponds to the case when ao = {U; 0; 0) 
and A{x) = in ([T|). We can also consider the situation when oq = and 
A{x) ^ in ([1]). In this case we note that if A{x) is a symmetric strictly 
positive matrix (e.g., A{x) = al, a > 0), then Lqv = A{x)v can be inter- 
preted as a drag/friction term which models the resistance offered by the 
fluid against its flow (see, e.g., p9] for some discussion). 

Thus, our general model includes the case of interaction of the Poiseuille 
(or Oseen type) flow (with a possible drag/friction) in the domain O bounded 
by the (solid) wall S and a horizontal boundary 17 on which a thin (non- 
linear) elastic plate is placed. The motion of the fluid is described by the 
3D Navier-Stokes equations linearized around the Poiseuille (or Oseen) flow 
oo(x). To describe deformations of the plate we consider a generalized plate 
model which accounts only for transversal displacements and covers a gen- 
eral large deflection Karman type model and can be also applied to nonlinear 
Berger and Kirchhoff plates (see the discussion in |14j and also in Section S]). 
Since we deal with linearized fluid equations the interaction model consid- 
ered assumes that large deflections of the plate produce small effect on the 
corresponding underlying flow. 

We note that the mathematical studies of the problem of fluid-structure 
interaction in the case of viscous fluids and elastic plates/bodies have a long 
history. We refer to [HI El El [201 EH |22l [231 [Ml EO] and the references therein 
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for the case of plates/membranes. The case of moving elastic bodies [T7] and 
the case of elastic bodies with the fixed interface [H HI |H [19] were studied; 
see also the literature cited in these references. We also mention the recent 
short survey [16j and the paper [15] which deals with dynamical issues for a 
model taking into account both transversal and longitudinal deformations. 
All these sources deals with the case of bounded reservoirs O. 

In this paper our main point of interest is well-posedness and long-time 
dynamics of solutions to the coupled problem in ([2])-([7]) for the velocity v 
and the displacement u in the case of unbounded domains O. 

In our argument we use the ideas and methods developed in our previous 
paper [14J. Our main difficulties in comparison with [14J is related to the 
facts that (i) we deal with the (possibly) unbounded domain O (hence, we 
loose some compactness properties of the fluid velocity variable and cannot 
use eigenfunctions of the Stokes operator) and (ii) the fluid equation ^ 
is perturbed by nonconservative and nondissipative term (hence, we can 
loose the energy monotonicity and need some additional argument for non- 
monotone parts). To overcome these difficulties we are enforced to use a 
general basis in the fluid component and a specially constructed extension 
operator Ext of functions on into solenoidal functions on O. 

The paper is organized as follows. In Section [2] we introduce Sobolev 
type spaces we need and provide with some results concerning the exten- 
sion operator Ext. In Section [3] we prove Theorem 13.21 on well-posedness in 
the case of linear model and study stability properties of solutions in The- 
orem 13.41 Section H] is devoted to nonlinear problem. We prove here that 
the problem generates a dynamical system (see Theorem 14. 3p which, under 
some additional conditions, possesses a compact finite-dimensional global 
attractor (Theorem 14. 6p . 

2 Preliminaries 

Let be a sufficiently smooth domain in and H^{D) be the Sobolev space 
of order s € M on D which we define (see [35]) as restriction (in the sense of 
distributions) of the space H'^i^'^) (introduced via Fourier transform). We 
define the norm in H^{D) by the relation 

= inf : w e H^W'-), w = u on L>| . 

We also use the notation || • ||d = || • \ \o,D for the corresponding L2 norm and, 
similarly, (•, ■)£, for the L2 inner product. We denote by Hq{D) the closure 
of Cq°{D) in H^{D) (with respect to || • \\s,d) and introduce the spaces 

Ht{D) := {/|^ : / G H'{R''), supp/ cd], s G M, 

to describe boundary traces on C dO. Since the extension by zero of 
elements from Hl{D) give us elements of H^{W^), these spaces Hl{D) can 
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be treated not only as functional classes defined on D (and contained in 
H^{D)) but also as (closed) subspaces of H'^{W^). We endow the classes 
Hl{D) with the induced norms = ||/||s,Rrf for / G H^{D). It is clear 

that 

< ll/li:,!,, f^Hl{D). 

It is known (see [351 Theorem 4.3.2/1]) that C^{D) is dense in Hl{D) and 
HI{D) = Wq{D) for -l/2<s<oo, s - 1/2 {0, 1, 2, . . .}. 

The norms || ■ ||^£) and || • \\s,D sre equivalent for these s. Note that in 

the notations of [27] the space H^^^^'^{D) is the same as H^^^^"^ (D) for 
every m = 0, 1, 2, . . . Below we also use the factor-spaces H^(D)/M. with the 
naturally induced norm. 

To describe fluid velocity fields we first introduce the class ^o(C') of C°° 
vector-valued solenoidal (i.e., divergence-free) functions v = (f^;f^;f^) on 
O which vanish in a neighborhood of dO and also for |x| large enough. Then 
we denote by X the closure of '^o(C') with respect to the L2-iiorm and by 
V the closure of '^o(C') with respect to the H^-norm. One can see that 

X = {v= {v^;v'^;v^) G [^2(0)]^ : divv = 0; jnV = iv,n) = on dO} 

and 

V CV"" = {v = {v^;v^;v^) G [H^{0)f : divv = 0; v = on dO} . (13) 

For some details concerning this type of spaces see, e.g., |25l|3l] and [18] . 

The following (geometry type) hypothesis plays an important role in our 
further considerations. 

Assumption 2.1 (Domain Hypothesis) We assume that 

(i) there exists a smooth bounded domain O' C O such that J7 C dO'; 

(ii) we have the equality in (fT3l) . i.e., V = V^. 

The sense of the first requirement in Assumption 12 . 1 1 is obvious. As for the 
second one we refer to [26] for a discussion of conditions on the domain which 
guarantee the equality V = (see also [THJ Sect.4.3] and the references 
therein). Here, as examples, we only note that this property holds in the 
following cases: (i) O is a smooth domain with the compact boundary; (ii) 
O = Mi = {x3 < 0}; (iii) O is given by ([9]) with smooth bounded ^ or 
with ^ as in (jl2p (infinitely long pipes and tubes of possibly varying cross 
section are also admissible). 

We also need the Sobolev spaces consisting of functions with zero average 
on the domain 0,, namely we consider the subspace 

L2{n) = |u G L2{n) : j u{x')dx' = o| 
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in L2{n) and also the subspaces H'{n) = H%n)nL2{n) in H'in) for s > 
with the standard H^{Q)-norm. The notations H^{il.) and i/Q($7) have a 
similar meaning. We denote by P the projection on Hq{Q) in Hq{Q,) which 
is orthogonal with respect to the inner product (A-, A-)q. As it was already 
mentioned in [14| the subspace {I—P)Hq{Q,) consists of functions u £ Hq{CI) 
such that A'^u = const and thus has dimension one. 

In further considerations we need the following assertion concerning ex- 
tension of functions defined on il,. 

Proposition 2.2 Let Assumption \2.lY i) he in force. Then there exists a 
linear hounded operator Ext : ^2(^1) 1— ?• [L2(C')]^ such that 

di-v Ext[ilj] = Q inO, {Ext[Tp],n)\^ = Tp, {Ext[ijj],n)\^ = 0, 

and 

ll^^*[V']ll[^v2-.(o)]' < cuu V5 > 0, y^e L2{n). 

Moreover, 

• ifil^e H'^{n) for some < s < 1, then Extiip] G [H'+^/^{0)f with 
the estimate 

and the relations Ext[ilj]\g = (0;0;0) and Ext['ifj]\^ = (0;0;'(/') on the 
houndary of dO; 

• there exists a smooth hounded suhdomain O' in O such that (i) Vt d 
dO', (ii) Ext[i}]\^^^, = 0, and (Hi) Ext[il:]\^, G [H'^{0')f' provided 

ij) G i?Q^^^'^(il) for some 5 > {). 

Proof. On a smooth bounded subdomain O' in O such that C dO' we 
consider the following Stokes problem: 

-v/\v + Vp = 0, div u = in 0'\ 
v = Q ondO'X^] w = (0;0;V) on f^, (15) 

where G L2{^) is given. This type of boundary value problems in hounded 
domains was studied by many authors (see, e.g., [25l[3l| and also the recent 
monograph |18j and the references therein). To construct an extension oper- 
ator we need the following properties of solutions to (fT5]) (for some discussion 
and references concerning the assertion below we refer to [14]). 

Proposition 2.3 Let ^ G Hl{Q) with -1/2 < s < 3/2 and J^'ip{x')dx' = 
0. Then prohlem (jl5p has a unique solution 

{v;p} G [H'+^/\0')f X [H'-^/^{0')/R] 

such that 

lbll[i?s+i/2(c)/)]3 + \\p\\h''-'>-/^{o')/r - collV'llHj(n)- 
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Now we can take a solution v to (jlSp and define Ext[ip\ as the zero extension 
of V on the domain O. One can see that for this operator Ext ah statements 
of Proposition 12.21 are in force. □ 

Remark 2.4 We could not find in the hterature an appropriate statement 
of Proposition 12.31 for unbounded domains. On the other hand we do not 
know an extension result in the class solenoidal functions with estimate 
()14p for some range of the parameter s. This is why we use this way for a 
construction of the operator Ext. We also note that in the case when O is 
bounded we can take O' = O. In this case Ext is a Green type operator 
which maps tp into v according to (llSp . Exactly this extension operator was 
used in [14j. 

Using the extension operator constructed above we introduce the spaces 
which we need to describe the interaction between fluid and plate. 
Let Assumption 12.11 be valid and 

M{0) = ^^v = VQ + Ext[il:] : G ^o(C), V' e -^o(f^)} ■ 

Then we denote by X the closure M{0) with respect to the L2-norm and 
by V the closure of M{0) with respect to the i?^-norm. One can see that 

X = {v = {v^]v'^]v^) G [^2(0)]^ : div?; = 0; 7„v = {v,n) = on 5} 

and 



V = VO = }v = {v']v';v') G [H\0)] 



div V = 0, u = on S*, 



We equip X with L2-type norm || • \\ci and denote by (•, ■)o the corresponding 
inner product. The space V is endowed with the standard norm. 

In conclusion of this section we mention that in the the case of the 
Poiseuille flow in the tube or between two planes described above we deal 
with a domain satisfying the Friedrichs-Poincare propertj{l|: 

3do>0: / \v{x)\^dx <d% f \Vv{x)\^dx, yveH^{0). (16) 
Jo Jo 



By the localization argument one can show that the inequality in (jlOp implies 
a similar property for any v £ {g G H^{0) : g\s = 0} and thus 

3co>0- \\v\\o < coW'^vWo, yveV, (17) 
for the Friedrichs-Poincare domains. 



^This property is valid in the case when the domain O is bounded at least in one 
direction. 



7 



3 Linear problem 



In this section we consider a linear version of ©-([HI) which is obtained by 
replacing equation with its linear counterpart. Thus we deal with the 
following problem 



vt — I'Av + Lqv + Vp = Gf{t) and divv = in Ox (0, +00) 
v = on S and v = {v^;v^;v^) = {0;0;ut) on Q, 
utt + A^u = Gpi{t) +p\n on il, 
du 



(18) 
(19) 
(20) 

(21) 



^_ -0 on / ut{x\t)dx' = d for all t > 0, 

which we supply with initial data of the form 

v{Q) = VQ, u(0) = uo, ut{G)=ui. (22) 

Similarly to [2] we consider the following class of test functions 

r <t^eL2[Q,T-[H\0)f), 0i GL2(0,r;[L2(O)]3), 
Ct=14> div</. = 0, ^\s = 0, <P\n = (0; 0; 6), 0(r) = 0, 
I 6 G L2(0,r;^2(5^))^ ^ L2(0,r;L2(J])). 

and introduce the following definition. 

Definition 3.1 A pair of functions {v{t);u{t)) is said to be a weak solution 
to the problem in (|18|) -(|22 p on a time interval [0, T] if 

. V G Loo(0,r;X)nL2(0,r;y); 

• u(^L^{Q,T-Hl{Q)), ut e Loo{0,T;L2{n)) and u{0) = uo; 

• for every (p € Ct the following equality holds: 



-/ {v,c^t)odt + v {Vv,V^)odt+ / {Lqv, 







Jo 







T 



{ut,bt)ndt+ I {Au,Ab)ndt 

T 



{Gf,<P)odt+ / (Gp,,6)Qdt + (7;o,</)(0))o + (ni,6(0))f^; (23) 

JO 

• the compatibility condition v{t)\fi = (0;0; nt(t)) holds for almost all t. 

The same argument as in [14J shows that a weak solution {v{t);u{t)) satisfies 
the relation 

{v{t),7p)o + {ut{t),/3)n = {vo,ip)o + {ui,f3)n 

v{Vv, Vip)o + {Lqv, tp)o 

+ (An, A/3)n - {Gf,i,)o - {GpuP)^] dr 



t r- 
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for almost all t € [0, T] and for all t/j = {tp^ ; ■ip'^ ; ip^) € W, where P = ''P^Iq 
and 

W = [ipev\4^\n = (0; 0; /3), /3 G ^^(f^) } • (24) 

It also follows from the compatibility condition and the standard trace the- 
orem that the plate velocity ut possesses an additional spatial regularity, 

1/2 

namely we have that ut G L2{0,T;H^' (Q)). 
Below as phase spaces we use 

n = ^{vq;uq;ui) e X X HQ{n) xL2{fl) : {vo,n) = = ui on (25) 



and 



-}{ = {{vo;uo;ui) e-H: uq e Hlin)^ CH (26) 



with the norm |j(uo; uq; 'Ui)||^ = \\vo\\o + II^'WoIIq + ll^illn- 

Our main result in this section is the following well-posedness theorem 
concerning the linear problem. 

Theorem 3.2 Let Assumption \2. 1\ be in force. Assume that 

Uo = {vo;uo;u^) en, Gf{t) e L2{0,T;V'), Gpi{t) e L2{0,T;H-'/\n)). 

Then for any interval [0,T] there exists a unique weak solution {v{t);u{t)) 
to ()18p - ()22p with the initial data Uq. This solution possesses the property 

U{t;Uo) = U{t) = {v{t)-u{t)-ut{t)) G C(0,T;X x Hl{n) x L2(J))), 

and satisfies the energy balance equality 

rt 



dT 



£o{v{t),u{t),ut{t))+ [ [iy\\Vv\\l + {Av,v)o 
Jo 

= £o{vo,uo,ui) + / {Gf,v)odT+ / {Gpi,Ur)ndT (27) 
Jo Jo 

for every t > 0, where the energy functional £q is defined by the relation 

£oiv{t),u{t),utit)) = \ {\\v{t)\\l + \\ut{t)\\l + \\^u{t)\\l) . 

Gf = and Gpi = 0, then Theorem 13.21 implies that the problem in (jlSp - 
(j22p generates a strongly continuous semigroup. In order to state our result 
on asymptotic stability of this semigroup we need additional assumptions. 

Assumption 3.3 (Stability Hypothesis) Assume that one of the follow- 
ing conditions is valid: 

• either the matrix ^4(2;) in ([1]) is uniformly strictly positive, i.e., 
3a >0: {A{x)C,Or^ > a\C\l,, G M^ x G O; 
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or the domain O satisfies the Friedrichs-Poincare property (|16p and 

3 ,5 > : {A{x)^, e)K3 > - - |eiM3 , Ve G M^ X G O, (28) 
where cq is the constant from the Friedrichs-Poincare inequality in 



Thus in the case of a general domain O satisfying Assumption 1 2 . 1 1 to obtain 
a result on long-time dynamics we need to assume the presence of some 
additional damping mechanism (drag/friction terms). If the domain satisfies 
the Friedrichs-Poincare property, then the result can be achieved without 
any damping (e.g., we can take A{x) = 0). Moreover, we note that the 
condition in (j28]) is true when sup^.^^ |A(x)| < vcq^, where |^(x)| is the 
operator (Euclidian) norm in M'^. In the case of the Poiseuille type flow (see 
(jlip ) this means that jVaj^^ajjaj is small enough. Since the profile a can be 
written in the form a = kv~^a^, where solves (jlOp with v = 1 and k = 1, 
the latter condition is satisfied when kv~'^ < c{B). Here k is the Poiseuille 
velocity parameter and c{B) is a constant depending on the cross-section B 
the tube O. In the case of the Oseen model we have gq = {U; 0; 0) in (fTTT) 
and thus there are no restrictions on the velocity U of the underlining flow 
for Friedrichs-Poincare domains. 

Theorem 3.4 In addition to the hypotheses of Theorem \3.S\ we assume that 
Assumption \3.3\ is in force. Then there exist positive constants M and'j such 
that for every initial data Uq = {vo;uo;ui) from % we have 

\\U{t)\\l < Me-^'\\Uofu + M re-^(*-) f||GKr)||^ + \\G,i{r)\ty,J dr. 

Jo ^ 

(29) 

In particular, if Gf = and Gpi = 0, then the Co-semigroup generated by 
()18p - (|22p is exponentially stable in H. 

In the case of a general operator Lq we need to add the term 

Mi + M2Sup\A{x)\] [ e-'^^*-^^\v{T)fdT 
xeo ) Jo 

in the right hand side of (|29p . Here 1^4(2;) | denotes the operator (Euclidian) 
norm in and Mi = when the domain O satisfies the Friedrichs-Poincare 
property in ([16]). 



Proof of Theorem lOl 

In the case when Lq = and O is bounded this theorem was proved in 
(see also [30] for a similar result). We use the same idea as in [14]. The main 
difficulty which we are faced is that we loose several compactness properties 
of the model (e.g., we cannot use the basis of eigenfunctions of the Stokes 
operator). 
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Step 1. Existence of an approximate solution. Let {V'jjigN be an (orthonor- 
mal) basis in the space V consisting of the smooth finite in O functions. 
Denote by {CijieN the basis in Hq{Q,) which consists of eigenfunctions of the 
following problem 

with the eigenvalues < ki < /t2 < . . . and \\£,i\\n = 1. Let = Ext[(,i], 
where the operator Ext is defined in Proposition 12.21 This proposition also 
yields (pi is in some vicinity of Q, and thus as in [H] one can conclude 
that (9x3 0f = on J7. 

We define an approximate solution as a pair of functions 

m n n 

1=1 j=l j=l 

satisfying the relations 

iVn,mit),x)o + {Un{t),h)n + Vn,'m{t) , ^x)o + ( Ati„ (t) , A/l)n 

= -{Lovn,m,x)o + {Gf{t),x)o + {Gpi{t),h)n (31) 
for t € [0, T] and for every x h of the form 

m' n' 

X = ^ Xki^k + Ext[h] with h = ^ hk^k, (32) 
k=l k=l 

where m' < m and n' < n. It is clear that x ^ ^ ^-iid x\q = (OjO; h). The 
system in ()3ip is endowed with the initial data 

Vv,m{0) = ^m{vo - Ext[ui]) + Ext[PnUi], 
Un{0) = PnPuQ + {I - P)uo, Un{0) = Pn^l, 

where 11^ is the orthoprojector on Lin{^j : j = 1, . . . , m, } in X and P„ is 
orthoprojector on Lin{^i : i = 1, . . . , n} in L2{^)- Since Ext : L2{^) i->- X, 
it is clear that 

{vv,mi0);un{0);un{0)) {vq;uq;ui) strongly in ?^ as m,n ^ oo. 

As in jH] one can show that (|3ip can be reduced to some ODE in M'"+"' 
and with given initial data has a unique solution on any time interval [0, T]. 
It follows from ([30]) that 

m 

(t) = ai{t)iPi + Ext[dtUn{t)]. 
1=1 

This implies the boundary compatibility condition: 

Vn,m{t) = {0;0;dtUn{t)) ou il. (33) 



11 



Step 2. Energy relation and a priori estimate for an approximate solution. 
Taking x = Vn^m and h = dtUn{t) in (j3ip we obtain the following energy 
balance relation for approximate solutions 

So{Vn,m{t),Un{t),dtUn{t)) (34) 
+ U / \VVn,m\'^dxdT + / {AVn,m, Vn,m)odT 

Jo Jo Jo 

= £o{Vn,m{0),Un{0),dtUn{0)) + {G f , Vn,m)odT + {Gpl, dtUn)Q.dT. 



We use here the structure of Lq which after simple calculations (see, e.g., 
Lemma 1.3 [Ml Ch.2]) yields the equality {LoVn,m,Vn,m)o = {Avn,m,Vn,m)o- 
The relation in (|34p and Gronwall's lemma implies the following a priori 
estimate 

sup {\\Vn,mmO + \\^Mt)\\l + WdtUnitml} 

ie[o,T] 

rT 



+ [ {\\VVn,m\\o + hn,m\\o)dT <Ct. (35) 

Jo 



By the trace theorem from ([33]) and ([35]) we also have that 
rT rT 

|2 



ll^t^in(r))||^y2^^^C?r = \\Vn,m{r)\\y^^Q^dT <Ct- (36) 

Step 3. Limit transition. By ([35]) the sequence {{vn,m', Un] dtUn)} contains a 
subsequence such that 

■,dtu n) {v]U]dtu) * -weakly in Loq(0,T;T-L)] (37) 

Un^u strongly in C{0,T; H^''{^})), Ve > 0; (38) 

Vn,m V weakly in L2{Q,T;V). (39) 

To obtain ([38]) we use the Aubin-Dubinsky theorem (see, e.g., [32l Corol- 
lary 4]). By ([36]) we can also suppose that 

dtUn dtu weakly in L2(0, T; hI'^{^1))] (40) 

Vn,m^v weakly in L2(0,T;/7i/2(aO)). (41) 



Applying the same argument as in |14j and using relations ([37p ~ ([4ip we con- 
clude the proof of the existence of weak solutions which satisfy the corre- 
sponding energy balance inequality. At this point we use Assumption 12. l( ii) 
to approximate elements from W by elements of the form (|32p . We need 
this to establish ([23]) for cj) G Ct- 

Step 4- Uniqueness. We first consider the case when Lq = and use 
Lions' idea (see ^28]), with the same test function as [14] in the case of a 
bounded domain. After establishing properties of solutions in this case we 
consider the term Lqv perturbation. 
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Let U^{t) = (v^ {t);u^t);ul{t)), j = 1,2, be two different solutions to 
the problem in question with the same initial data and Lq = 0. Then their 
difference U{t) = U^{t) — U'^{t) = {v{t);u{t);ut{t)) satisfies the variational 
equahty 

T /-T f-T nT 

{v,(t>t)o + v {Vv,V(l))o- {utM)n+ / (An,A6)Q = 

JO JQ 

for all (f) G Ct, b = Now for every < s < T we take 

{ - I dr [ dCv{C), t < s, 
^ ^^(t) = 1 Jt Jo 

(o, t>s, 

as a test function. The same calculation as in jl4j yields the uniqueness in 
the case Lq = 0. 

Step 5. Continuity with respect to t and the energy equality. Using 
the Lions lemma (see [271 Lemma 8.1]) by the same argument as in [14J 
we first prove any weak solution {v{t);u{t);ut{t)) is weakly continuous in 
X X Hl{n) X L2(0). 

To prove the energy equality (in the case Lq = 0), we follow the scheme 
of [28l Ch.l], see also [27, Ch.3], in the form presented in [H]. Thus as in 
|14j we can conclude that the solution is strongly continuous in t. Moreover, 
the energy relation in the case Lq = with Gj = and Gpi = implies that 
the corresponding solutions generates strongly continuous semigroup. 

Step 6. Case Lq ^ 0. Using the energy relation for the problem with 
1/0 = and Gf{t) := Gf(t) — LQv{t) we can establish the uniqueness of solu- 
tions via the Gronwall's type argument and also the smoothness properties 
in the general case. This completes the proof of Theorem 13.21 



Proof of Theorem 13.41 

To prove the estimate in ()29p . we construct a Lyapunov function using an 
idea from [8J (see also p^). Let 

V{vo, uo, ui) = £o{vo,uo, ui) + e^'(?;o, -uo, ""i), 

where ^'(fo, uo,ui) = (uq, + {vq, Ext[uo])o and e > is a small param- 
eter which will be chosen later. We consider these functionals on approxi- 
mate solutions {vn,m',Un) for which Puq = uq and thus Pun{t) = Un{t) for 
all t > 0. This allow us to substitute in ()3ip Ext[un] instead of x and obtain 
that 

d d 

^*n,m(t) = -^'^{Vn,m{t),Un{t),dtUn{t)) 

= \\dtUn\\h Ext[dtUn])o - {LoVn,m-,Ext[Un])o 

- h'{VVn,jn,^Ext[Un])o " ll^^^nUn 

+ {Gf,Ext[Un])o + {GpUUn)n- (42) 
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By Proposition 12.21 using the compatibility condition in and the trace 
theorem we have that 

\{Vn,m,Ext[dtUn])o\ < C\\Vn,m\\o\\dtU„\\n < C [\\VVn,m\\o + lkn,m||c)] • 

Similarly, for every r/ > we have 



VExt[Un])o\ < 'nW^UnWn + Cr, [\\VVn 



m\\0\ 



and 



G/lly/ + \\G -"^ 



\{Gf,Ext[un])o + iGpi,Un)n\ < r?||An„||^ + 
It is also clear that 

\{LoVn,m, Ext[Un])o\ < ??ll^'"n||n + Crj [II Vl'n^mllc) + ||^n,r>i||o 

that 



Therefore it follows from 

d 
di 



112 1 



+ c 



GfWyl + IIG' 



pi\\-i/2,n 



Using the energy relation in ()34p we also have that 

+ ^||'*^n,m||c) 



GfWyl + 116*' 



pi\\-i/2,n 



Vn,m)o, Vr/ > 0. 



One can see that the function Vn,m(i) = y{vn,mit),Un{t),dtUn{t)) satisfies 
the relations 

ao£o{Vn,m{t),Un{t),dtUn{t)) < K,m(i) < ai£o{Vn,m{t) , Unit) , dtUn{t)) 

for sufficiently small e > 0. Using the stability hypothesis in Assumption 13. 
we can choose rj > and o" > such that 



{l^ - r])\\VVn,m\\o - 'n\\Vn,m\\o + i^'"' 

Therefore we have that 
d 



n,m7 ^n,m 



> CT [\\VVn 

,ni llo 

+ 



dt 



Vn,m{t) + a2Vn,m{t) < 03 



I Gf 1 1 y / + 1 1 G. " ^ 



pi\\-i/2,n 



with positive constants Oj. This implies relation (|29|) for approximate solu- 
tions. The limit transition yields (|29p for every weak solution. 

In the general case we can apply (p9]) with Lq := (oq, V)v+hv and Gf := 
Gf — Av + fiv, where /i > (in the case of the Friedrichs-Poincare domains 
we can take = 0). This implies the desired conclusion and completes the 
proof of Theorem 13. 4[ 
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4 Nonlinear problem 

In this section we deal with problem (I2])-® with a nonlinear feedback force. 
First we impose the following hypotheses concerning the force J-{u) in the 
plate equation 

Assumption 4.1 • There exists e > such that J-{u) is locally Lips- 
chitz from Hq^'^^Q) into in the sense that 



for any u-i € Hq{^) such that ||uj||2,n ^ ^• 

• There exists a -functional n(ti) on Hq{Q,) such that J-{u) = n'(u), 
where H' denotes the Frechet derivative ofH. 

• The plate force potential 11 is hounded on hounded sets from HQ{il) 
and there exist r] < 1/2 and C > such that 



Examples of nonlinear feedback (elastic) forces J-'{u) satisfying Assump- 
tion HTT] are described in [9j and [Mj, see also [13 1 . They represent different 
plate models and include Kirchhoff, von Karman, and Berger models. 

4.1 Well-Possedness 

Definition 4.2 A pair of functions {v{t);u{t)) is said to be a weak solution 
to ©-(HD on a time interval [0, T] if 

. v£L^{0,T;X)f]L2{0,T;Vy, 

• ue L^{0,T;H^{n)), ute L^{0,T;L2{^)), u{0) = uo; 

• the equality in 1^3^ holds with Gpi{t) := -T{u{t)) + Gpi(t); 

• the compatibility condition f(t)|n = (0;0; nt(t)) holds for almost all t. 
Theorem 4.3 Let Assumptions \2. i] and \4-l\ he in force. Assume that Uq = 



{vo;uo;ui) e H, Gf{t) G L2{0,T;V') and Gpi{t) G L2{0,T; H'^/^in)). 



Then for any interval [0,T] there exists a unique weak solution {v{t);u{t)) 
to ([2]) -([8]) with the initial data Uq. This solution possesses the property 



where T-L is given hy [25\) . and satisfies the energy halance equality 



\\T{ui) - J^2{U2)\\^1/2,Q < Cr\\ui - U2\\2-e,n 



(43) 



r/||An||^ + n(u) + C > , G /7^(0). 



(44) 



U{t) = {v{t);u{t);ut{t)) £ C{0,T;n) 



(45) 



£{v{t),u{t),ut{t))+ / [u\\Vv\\l + {Av,v)o]dT 




£{vo,uo,ui) + / {Gf,v)odT+ / {Gpi,Ur)ndT (46) 
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for every t > 0, where the energy functional E is defined by the relation 
£{v,u,ut) = ^ {Wvfo + htWfi + ||An||^) +n(u). 

Moreover, there exists a constant ap^ x > such that for any couple of weak 
solutions U(t) = {v{t);u{t);ut{t)) and U{t) = {v(t);u{t);ut{t)) with the 
initial data possessing the property ||?7o||'H7 l|t^o||'H < R we have 

r||V(^^-*)llo^^<ai?,T||f/o-^7o|||^, tG[0,r]. (47) 
Jo 

The spatial average of u{t) is preserved. In particular, if Uq G 7i, then 
U(t) for every t > 0. We recall that % is defined by [2S\) . 

Proof. The proof of the local existence of an approximate solution is al- 
most the same, as in the linear case (see Theorem I3.2p . We use approximate 
solutions of the same structure which satisfy ([3T]) with —F{un{t)) + Gpi{t) 
instead of Gpi(t). Then using the standard argument we establish the energy 
relation in ()46p for these approximate solutions. Now the positivity type es- 
timate in (I44p allow us to obtain the same a priori estimates as in (I35p and 
(j36p . Therefore we can prove the global existence of approximate solutions 
and establish the existence of a weak solution U(t) = {v{t);u{t);ut{t)) by 
the same argument as in the linear case. To make limit transition in the 
nonlinear term we use (j43p . 

Next we consider the pair {v(t);u{t)) as a solution to the linear problem 
with Gpi{t) := —T{u{t)) + Gpi{t). This allow us to obtain (j45p and also 
derive energy balance relation (146p from (127p using the potential structure 
of the force T: J^{u) = Yi'{u). 

Since the difference of two weak solution can be treated as a solution to 
the linear problem with Gf = and Gpi{t) := T{u{t)) — F{u{t)), we can 
obtain (I47p from the energy equality (I27p . The uniqueness follows from (I47p . 

Preservation of the spatial average of u{t) follows from the same property 
for approximate solutions. □ 

We can derive from Theorem 14.31 the following assertion. 

Corollary 4.4 In addition to the hypotheses of Theorem \4-3\ we assume 
that Gf{t) = Gq € V is independent oft and Gpi{t) = 0. Then problem 
©^([H]) generates dynamical systems (StjIL) and {St^H) with the evolution 
operator defined by the formula SiUq = {v{t);u{t);ut{t)), where {v;u) is a 
weak solution to ({2]) -([8]) with the initial data Uq = {vo;uo;ui). 

If we assume in addition that Go = and Assumption \3.3\ holds, then 
these systems are gradient with the full energy £{vo,uo,ui) as a Lyapunov 
function. This means that (a) U i— )• £{U) is continuous on %, (b) £{StUQ) 
is not increasing in t, and (c) if £{StUo) = £{Uq) for some t > 0, then Uq 
is a stationary point of St (i.e., StUo = Uq for all t > 0). Moreover, the set 
£r = {Uo ■ £{Uq) < R} is a bounded closed forward invariant set for every 
R>0. 
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Proof. The argument is the same as in [14) . We only note that under 
Assumption 13.31 from (j46p (with Gf = and Gpi = 0) we have that every 
stationary point C/* for St has the form [/* = (0; u; 0), where u G Hq{Q). □ 

4.2 Stationary solutions 

As above we assume that Gpi = and Gf = and Assumptions 12.11 and 13.31 
holds. It follows from Definition 14.21 that a stationary (time-independent) 
solution is a pair {v; u) from V x Hq{Q,) satisfying the relation 

i^iVv, VV')o + (^0^^, ^)o + {Au, A/3)n + iHu),/3)n = (48) 

for any Tp £ W with ^'^|q = where W is given by ()24p . Taking ip = v 
we conclude that i/||Vf ||^ + {Av,v)c) = and hence from Assumption 13.31 
we have v = Therefore we obtain the following variational problem for 
u € Hl{n): 

(A^x,A/3)o + (^(^x),/3)f, =0, V/3G^o'(f^)- (49) 

As in [14] we can show the existence of a family of solutions to (I49p param- 
eterized by a real parameter. In the case of the zero average of u we can fix 
this parameter and obtain the following assertion (see [14] for details). 

Proposition 4.5 ( |14| ) In addition to Assumption we assume that 
there exist r] < 1/2 and c > such that 

T]\\Au\\l + {u,T{u))n>-c, yu£H^{n). (50) 

Then the set TVq of solutions u to problem with the property udx = 
is nonempty compact set in Hq{Q). This implies that the set J\f of all 
stationary points of St in the space % is nonempty compact set and has the 
form 

A/'= |(0;'u;0) : -u € -&^(J1) so/ves g^j . (51) 

4.3 Asymptotical behavior 

In this section we are interested in global asymptotic behavior of the dy- 
namical system {StjH). Our main result states the existence of a compact 
global attractor of finite fractal dimension. 

We recall (see, e.g., [3l[71|33]) that a global attractor of the dynamical 
system {St-,'H) is defined as a bounded closed set 21 C 'H which is invariant 
{S{t)% = 21 for all t > 0) and uniformly attracts all other bounded sets: 

lim sup{dist-^(S'(t)y, 2t) : y € B} = for any bounded set B in T-L. 

Theorem 4.6 Let Assumptions \2.I[ \3.3l and \4-l\ be in force. Assume that 
Gpi = 0, Gf = and ^5^) hold. Then the dynamical system {St^T-L) possesses 
a compact global attractor 21 of finite fractal dimensioi^. 
Moreover, 

^ For the definition and some properties of the fractal dimension, see, e.g., [7] or [33] . 
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(1) Any trajectory 7 = {(v{t);u{t);ut{t)) : t G M} from the attractor 21 
possesses the properties 

{vf, uf, utt) G Loo(M; X X Hl{n) x L2{n)) (52) 

and there is R > such that 

supsup{\\vt\\o + \\ut\\l^ + \\uu\\n) <R^- (53) 

(2) The global attractor'^ consists of full trajectories {{v{t);u{t);ut{t)) : 
t € M} which are homoclinic to the set N , i. e. 

lim inf {\\v{t)\\% + \\u - u4l ^ + \\ut\\li) = 0, 
where A/q = |u G ^^(0) solves ig)}. /n addition we have 



lim dist^(5iy, AA) = for any initial data y ^H. (54) 

We emphasize that Theorem 14.61 deals with dynamics in the space T-L (the 
case of the zero spatial average of the deflection). For a possible approach 
to description of the system long-time behavior in the space T-L we refer to 
[Til Remark 4.9]. 

To obtain the result stated in Theorem l4.6l it is sufficient to show that the 
system is quasi-stable in the sense of Definition 7.9.2 [H] (see also Section 
4.4 in |12j). For this we can use the stability properties of linear problem 
(jl8p - (|22p established in Theorem [331 and the same argument as in [H] which 
yields the following assertion. 

Lemma 4.7 (Quasi-stability) Let the hypotheses of Theorem \4-6\ be in 
force and U''-{t) = {v^{t);u'^{t)]u\{t)), i = 1,2, be two weak solutions with 
initial data Uq = {vq\Uq\u\) from % such that \\Uq\\-h < R, i = 1,2. Then 
the difference 

z{t) = u\t) - u\t) = {v{ty,u{ty,ut{t)) 

satisfies the relation 

\\Z{t)\\l^<Mne-^''\\Zofn + Mn /* ||n(r)||2dr (55) 

Jo 

for some positive constant Mr and 7*. 

Proof. See fll] for some details. □ 
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To complete the proof of Theorem 14. 61 we note that by Proposition 7.9.4 [llj 
{Stil-L) is asymptotically smooth, i.e., for any bounded set B mT-L such that 
StB C B for t > there exists a compact set K in the closure B oi B, such 
that SfB converges uniformly to K. By Corollary 14. 41 the system is gradient. 
Proposition 14.51 yields that the set Af of the stationary points (see ([5T]) ) is 
bounded in 7i. Therefore to prove the existence of a global attractor we can 
use well-known criteria for gradient systems (see, e.g., [31 , Theorem 4.6] or 
Corollary 2.29 in [10]). 

The standard results on gradient systems with compact attractors (see, 
e.g., O [71 [33]) imply ([51]) . Since {St^T-L) is quasi-stable, the finiteness of 
fractal dimension dimjSl follows from Theorem 7.9.6 |11] . To obtain the 
result on regularity stated in (|52]) and (f53]l we apply Theorem 7.9.8 [TT] . 
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